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Abstract 

The rapid deposition of energy by Edge Localised Modes (ELMs) onto plasma facing components, 
is a potentially serious issue for large Tokamaks such as ITER and DEMO. The trigger for ELMs is 
believed to be the ideal Magnetohydrodynamic Peeling-Ballooning instability, but recent numerical 
calculations have suggested that a plasma equilibrium with an X-point - as is found in all ITER-like 
Tokamaks, is stable to the Peeling mode. This contrasts with analytical calculations (G. Laval, R. 
Pellat, J. S. Soule, Phys Fluids, 17, 835, (1974)), that found the Peeling mode to be unstable in 
cylindrical plasmas with arbitrary cross-sectional shape. However the analytical calculation only 
applies to a Tokamak plasma in a cylindrical approximation. Here, we re-examine the assumptions 
made in cylindrical geometry calculations, and generalise the calculation to an arbitrary Tokamak 
geometry at marginal stability. The resulting equations solely describe the Peeling mode, and 
are not complicated by coupling to the ballooning mode, for example. We find that stability is 
determined by the value of a single parameter A' that is the poloidal average of the normalised 
jump in the radial derivative of the perturbed magnetic field's normal component. We also find 
that near a separatrix it is possible for the energy principle's 6W to be negative (that is usually 
taken to indicate that the mode is unstable, as in the cylindrical theory), but the growth rate to 
be arbitrarily small. 

PACS numbers: 52.55.Tn,52.30.Cv,52.55.Fa,52.35.Py 



* Electronic address: |anthony. webster@ukaea.org.uk 



1 



I. INTRODUCTION 



Thermonuclear fusion requires plasmas with a pressure of at least an atmosphere, and 
temperatures in excess of 100 million degrees Kelvin. These conditions can be achieved 
in Tokamaks such as JEt[i], but the plasmas are subject to a number of instabilities, the 
consequences of which range from benign to structurally damaging. By understanding the 
instabilities that can occur, they can be avoided or mitigated. A class of instabilities that 
are only partly understood are Edge Localised Modes (ELMs) [2] . ELMs can lead to a rapid 
deposition of energy onto plasma facing components, and this is a potentially serious issue 
for proposed large tokamak devices such as ITErQ]. 

Our present understanding of ELMs is based on the linear ideal Magnetohydrodynamic 
Peeling-Ballooning instability (Wilson et al [4], Gimblett et alj^), which is thought to trigger 
ELMs, that subsequently evolve non-linearly. The studies upon which this understanding 
were based considered Tokamak equilibria with a smoothly shaped magnetic flux-surface at 
the plasma- vacuum boundary. In contrast, modern Tokamak plasmas have a cross-section 
in which the outermost flux surface is redirected onto divertor plates, forming a separatrix 
with a sharp "X-point" where the magnetic topology changes from closed (confined plasma), 
to open field lines along which plasma can flow to the divertor plates l|] . 

The first numerical evidence for a stabilising effect from the separatrix was found by 
Medvedev et al[6j. More recently, numerical studies of the Peeling-Ballooning instability 
in these X-point plasmas (HuysmansQ), have found that as the plasma's outermost flux 
surface is made increasingly close to that of a separatrix with an X-point, the Peeling mode 
becomes stabilised. Crucially the stabilisation appeared to happen before the plasma formed 
a separatrix with an X-point, shaping alone appeared to be sufficient to stabilise the mode. 
This appears to be contrary to theoretical work by Laval et al[8!], which has indicated that 
the peeling mode is unstable in cylindrical plasmas with an arbitrarily shaped cross-section. 
In addition the ELITE code (Wilson et alfl) has recently been used to examine Peeling mode 
stability as the outermost flux surface approaches the separatrix. It was found that although 
the growth rate reduced in size as the boundary more closely approximated a separatrix, it 
did so increasingly slowly, and its asymptotic behaviour was uncertain (Saarelma[l^). To 
help understand and reconcile these results, here in the first part of this two part paper we 
re-examine the assumptions made in the derivation of the peeling mode stability criterion 
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for a cylinder, and generalise the calculation so that it applies to a toroidal Tokamak plasma. 

As with the original studies of the Peeling mode in a cylinder, for simplicity we will 
firstly consider marginal stability, and generalise the condition for Peeling mode stability 
in a straight cylinder to a condition for Peeling mode stability in an arbitrary cross-section 
Tokamak plasma. The resulting equations only describe Peeling mode stability, and are not 
complicated by coupling to the Ballooning mode instability, for example. 

This paper generalises previous analytic calculations in a number of ways. Firstly it 
applies to axisymmetric toroidal geometries, as opposed to the cylindrical geornetry in w lich 
the Peeling mode has been extensively studied (for example see Laval et al[8], Lortz isl] . 
Connor et al[14!|). It allows for equilibrium poloidal currents at the plasma edge, in addition 
to the toroidal current that is solely included in previous analytic studies. The skin currents 
that are induced by a plasma perturbation are related to the difference between the magnetic 
field in the plasma and the vacuum, and it is found that at marginal stability a plasma 
perturbation induces a skin current that is parallel and proportional to the equilibrium 
edge current and proportional to the amplitude of the radial plasma displacement. The 
complicated-looking plasma-vacuum boundary condition that is usually found in association 
with the energy principle may be expressed as a simple relationship between the normal 
components of the plasma and vacuum magnetic fields. This is used to relate the generalised 
equations for Peeling mode stability at marginal stability to the energy principle, and this 
allows us to define the Peeling mode in terms of the energy principle's 6W. With this 
energy principle for the Peeling mode we can consider the trial function used by Laval et 



al 



finding that a single parameter A' determines the sign and magnitude of 6W. Finally 



the instability's growth rate is considered. 



II. BACKGROUND 

Laval et aljs^ considered a large aspect ratio ordering that neglects toroidal effects, and 
also neglects any equilibrium poloidal current at the plasma edge. The work suggests that 
the Peeling mode will be unstable for a non-zero edge current, regardless of the plasma cross 
section. Later we will reconsider Peeling mode stability for arbitrary cross section Tokamak 
plasmas, but firstly we consider some properties of the trial function considered by Laval et 
al 
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FIG. 1: A plot of the trial function used by Laval et al,8,], that for the most unstable modes 



with m nq has ~ ginqe ^ with 9 the usual straight field line angle 



In the second part to 



this paper we will show that the length / along a flux surface in the poloidal cross section may be 
parameterised by a, with a = — 7r..7r, as can 9. It is also shown that a physically reasonable model 



for 9 has 9{a) ^ J and /(a)//(vr) ^ (^f + 



1/4 



1/4 



da. Therefore 



daj § 

by parameterising both 9{q) and /(a), then the figure plots 9 (vertical axis) versus 1/1{tt) (horizontal 
axis), for e = 0.001 and n = 20. 



Note that for an element of length dl along a flux surface in the poloidal plane, 
= ^^^dx = ^^dx with the Jacobian, Bp the poloidal field, u = ^ is the 
local field-line pitch, and 0? orthogonal toroidal co-ordinate system (for ex- 

ample, see Freidberg[l2| for details). Thus the poloidal angle used in Laval et al, 
27r J^/ ^ = 27r J^dx' / § J^dx' = 6* is the same as the usual straight field line angle 11 1, 
and Q = ^ § ^dx = § The perturbation they consider has a plasma displacement 

^ ~ e*"*^, so if we plot ^ versus the length along a flux surface in the poloidal plane, then near 
the separatrix in an X-point equilibrium ^ will oscillate arbitrarily rapidly as we approach 
the X-point (see figured]). The most unstable modes have m ~ ng, so in the figure we plot a 
mode with m ^ rag, for which e*"*^ ~ ^mj^vdx^ Alternately, when m <^ nq or m ^ nq, then 
the mode localises poloidally in the vicinity of the X-point, and is approximately constant 
elsewhere. The rapid oscillation of ^ near the X-point makes it questionable whether it is 
physically acceptable, and other terms beyond ideal MHD need to be considered, but it cer- 



tainly means that the closer we approach the separatrix the greater the number of Fourier 
modes required (since m ~ nq), and the smaller a computer code's mesh spacing would need 
to be to represent the mode. Therefore as we approach the separatrix it will be increasingly 
difficult for a numerical calculation to represent the mode. 



III. CYLINDRICAL PLASMAS 



Here we outline the derivation of the marginal stability condition for an arbitrarily large 
aspect ratio (cylindrical) equilibrium with the Tokamak ordering (Freidberg 12j, Wesson jj). 
We start from the usual force balance equation J AB = Vp, and take the curl of both sides, 
expanding to give, 

= B.VJ- J.VB (1) 
Linearising the equation then gives for the equilibrium quantities 



= Bo.VJo - Jo.Vfio 



(2) 



and for the perturbed quantities 



= Bo.WJi + 5i.VJo - /o-V5i - Ji.V^o + (fio-V/o - Jq.VBo 



(3) 



Because Eq. |2] holds everywhere, with Bq.VJq — Jq.VBq having the constant value of zero, 
the last term in Eq. [3] that arises from the displacement of the plasma surface by ^, is zero. 
In the large aspect ratio approximation (Freidberg Wesson 1]), Eq. [^further simplifies 

to 

= 5o.VJi + 5i.VJo (4) 

Again, in the large aspect ratio (cylindrical) approximation we may write -Bi = A VV', 
for which 

= V A (^e; A VV^) 



(5) 



With this same ordering (Freidberg[l2|, Wesson[l|), Jq is taken to be parallel to e^. Therefore 
the Cz component of Eq. H] gives 



= BqV Ji.e, + 5i.VJ( 



'Oz 



(6) 



5 



Substituting for Ji then gives 



= 5o.V V'^A +5i.e 



dJ( 



Oz 



dr 



(7) 



Symmetry with respect to both the axial and the poloidal coordinates, means that we 
need only consider a single mode, 



ikz+imd 



(8) 



where r, 6', and z are the usual cylindrical coordinates, is a dimensional mode number, 
and m is a non-dimensional poloidal mode number. Hence we have 



ikz+imd 



im 

ikcz H 1 

r 



dr 



and because Bi = A Vtp, then Bi.Vr = —eg.Vtp, giving 



r 



We also have that 



m 



1 d dipn 



ikz+imd 



r dr dr 



Hence Eq. [7] gives 

= ( ikB, + im^ 



Br 



m 



1 d dipr^ 



-k 5- W'm + -TT^ 



r dr dr 



dJo4 
dr 



(9) 



(10) 



(11) 



(12) 



V.-B = requires the normal component of B to be continuous across the plasma-vacuum 
surface, that with Bi = Cz A Vip requires ipm to be continuous across the surface. The 
perturbed field from dipm/dr may be discontinuous however. Also Jqz = Ja just inside the 
plasma, but Jqz = in the vacuum outside the plasma, so if we integrate from an arbitrarily 
small distance inside the plasma surface to an arbitrarily small distance outside the surface, 
we get 







Bp 

ikBz + im — 
r 



dr 



1pm [\Jo. 

r 



(13) 



where [|/|] indicates the difference between / evaluated in the vacuum just outside of the 
plasma, and / evaluated just inside the plasma ([|/|] does not equal the integral of / from 
just inside to just outside the plasma). Writing k = — ^ and q = rB^/RBp, then gives 







m — nq 



m 



, dlpm 

dr 



+ 



B„ 



(14) 



as the condition for marginal stability. 



IV. TOKAMAK PLASMAS 



A. Assumptions 

Here we re-examine the assumptions made in going from Eq. [T2]to Eq. [T31 To obtain Eq. 
we allowed the perturbed fields to be discontinuous, but required the equilibrium fields 
to be continuous across the plasma- vacuum boundary. Because a discontinuous magnetic 



field requires a skin current [l3l]. we will allow perturbed skin currents, but the continuous 
equilibrium fields imply zero equilibrium skin currents. Therefore we assume that: (a) 
there are no equilibrium skin currents, but Jo can be discontinuous at the plasma-vacuum 
interface, and (b) perturbations to the magnetic field induce surface skin currents. The first 
of these assumptions means that with the exception of Jo, the equilibrium quantities will be 
continuous across the plasma-vacuum boundary. So if we integrate along the unit normal 
to the plasma surface, from a distance e just inside the surface to a distance e just outside 
the plasma, then we will get 

f{7i{l'))dl' = tf{ii{l)) ^ as e ^ (15) 

Ji-t 

where Tjil') is a path parameterised by /', that is parallel to the unit normal to the plasma, 
and where / is a continuous equlibrium quantity. We also get 

i ^^^' = /(^ + e)-/(/-e)^Oase^O (16) 

However, because the equlibrium current is discontinuous at the plasma surface being non- 
zero within the surface but zero in the vacuum, its derivatives act like delta functions, and 
hence 

fdr = Jo(/ + e)-Jo(l-e)= Jo (17) 



Jl-e 

and similarly 



f^dl' = fil)[\jo\] (18) 

The second of these assumptions, that the perturbation will produce skin currents at the 
plasma surface (i.e. that the perturbed magnetic field can be discontinuous at the perturbed 
plasma surface), means that for a perturbed current Ji 



i.e. the skin current acts like a delta function at the plasma surface. The reason for including 
the factor of 1/ RBp will be clear later. Similarly if / is a continuous function at the plasma 
surface, 

r fUl')dl' = f{l)^ (20) 
A few other remarks are worthwhile. Firstly the terms in Eq. [3] are evaluated at the 
equilibrium surface positions, but give the value of the perturbed force balance equation at 
the perturbed surface. Secondly, the perturbed unit normal fi = fio + fii, with ni ~ |^| <^ 1. 
Hence we have that 

— Jl-e J R2B^ 

= //-;7t(^-v/o-ni.vjo) di' 
= i!-:^n.vjo+o{o 

Jl-e RB„ dU 



(21) 



RB„ 



or equivalently. 



where J^^ indicates an integral from just inside the last closed flux surface, to just outside 
it. Similarly because Vip.n = Vip.nQ + 0{C,) = RBp + 0{^), then 

'^'^^^ = 1^ r fd^ + OiO (23) 
and 

l"ip+ 

a= / Ji# + 0(e') (24) 

Bearing the above remarks in mind, we now integrate Eq. [3] across the plasma surface, 
distinguishing perturbed quantities by a subscript of 1 and the equilibrium currents and 
magnetic field by a subscript of 0, to get 

= Bo.V j^J Jidij + j^^ (Si.V^) - (^j^J J.dij^ .VBo + O (e) (25) 

Note that because the terms in Eq. [3] are of order ^ (it describes a linearised perturbation 
to the plasma), the O(^) corrections that arise when integrating along the normal to the 
surface produce terms of order and are neglected. Therefore at leading order we have 



= Sq.Vcti + Bf 



ai.V^o (26) 
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where Bf = Bi.Vip. This is the generahsed force balance equation for Peehng mode marginal 
stability, valid for an arbitrary cross-section Tokamak plasma. The last term is a new term 
that is not present in a circular cross-section cylindrical geometry. 

The procedure of integrating across the plasma- vacuum boundary is clearer if we project 
out the components before integrating across the surface. We have done this as a check, but 
it is algebraically cumbersome, and obscures the physical arguments that are clearer in the 
presentation above. 



B. Ampere's law at the surface 



Before proceeding it is worth examining some relationships between the magnetic fields 
and the skin currents. At the plasma-vacuum interface, V.-B = and Ampere's law imply 



that 



13| 



n.B 

nA{Bv-B 



a 



RBr, 



(27) 
(28) 



with a = J^^ Jidip as before, n denotes the unit normal to the surface, and By the magnetic 
field in the vacuum. Because we assume zero equilibrium skin currents 



RBr, 



noA{B^-B, 







(29) 



and therefore because Eq. [57] gives Uq. yB^ — Bqj = 0, it then follows that Bq — Bq is not 
parallel to no, and hence the only way to satisfy Eq. [29] is if Bq = Bq. 
Considering the lowest order perturbation, we have 



(30) 



^ =niA[B^ -Bo] +noA[B^ -B, 



+eV no A [B^ -fi, 



which using Bq = Bq and Vip = RBpfio, gives 



a = ViIj A [BY - Bi 



and hence 



V^.a = 



(31) 



(32) 
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In the ijj, Xi coordinate system it is possible to directly evaluate a = //j*^^ V A Bidl', 
giving 

a = R^V(j)[\R,.B,\] - R^Bp[\V(l)A\] (33) 



C. Skin currents at marginal stability 



Using Vip.a = 0, while projecting out components of Eq. 

= Bq.V (a.Bp^ + Bi.Vi) Bp.Jo 
= Bq.V (a.V0) + Bi.Wip \ V0. Jo 



gives 



a.VB'^ 



Because 2J^ = -a.V {-j^) 



+ 2/V 

'^-^^ Bp.V (;^), then we can rewrite Eq. 







= 5o.V (a.V0) + 5i.VV^ 



V0.Jo 



Next we note that for Bq = + V0 A Vip, then 



givmg 



V0.Jo 



Bp.Jo 

V0.Jo 



Va 



H2 



Using Bi.Vip = Bq.V^^ and also noting that / = I{ip), Eqs. [35] now become 








5o. V ( ^ + 



(34) 



as 



(35) 



(36) 



(37) 



(38) 



Hence we have that 

^ = +f{4>-r ^dx') 

where / is some function oi (p — ^dx', so that Bq.V f = 0. Therefore using Eq. 
write Eqs. |39]in terms of the equilibrium current Jq, we have the skin currents at mar, 
stability given by 

^ = Jo^i, + Bof ((p- [ J^dx' 



(39) 



to 
inal 



(40) 
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— * 

If we require that = for = 0, then because Bq.V^^ 7^ 0, we simply have that 

a = Joi^ (41) 

i.e. that the skin current due to a perturbation at marginal stability is equal to the product 
of the equilibrium current at the edge and the radial displacement of the plasma. 

Notice that at this point we have not used any expansion, e.g. in terms of straight field 
line co-ordinates, and we have not neglected any poloidal dependencies of the equihbrium 
quantities. We have implicitly assumed the existence of the 0, ip coordinate system, but 
it is possible to re-express the coordinates in terms of arc length along a flux surface, and 
the consequent relations expressed in Eqs. [33] for example, remain valid (except at the point 
of zero size that is the X-point). Hence the results appear valid at the separatrix. 



D. Marginal stability 

To relate this to the previous cylindrical condition for marginal stability, Eq. [TU we 
consider V.-B, which may be written 



1 / d 



d 



d 



(42) 



with the Jacobian. Then we evaluate [|V.-B|], the difference between V.-B evaluated in 
the vacuum just outside the plasma and V.-B evaluated just inside the plasma. 



Noting that J^Vx-Bi = J^j^^^.Bi = then [|V.-Bi = 0|] requires that 



1 d 

— — J^V^.B, 



B2 } 

V 



Jx^X 



Bp.Bi 



52 

p 



d 



V0.-Bi 



(43) 



Eq. 



gives 



_ a. Bp 

a.V0 = 



V0.-Bi 
Bp.Bi 



(44) 



So we have 







Jxdx V B 



d 



Br, 



R^Bl 



(45) 



To simplify this we consider a limit of high toroidal mode number n and note that 
j-^ = -Bq.V — Then we note that -Bq.V is of order 1 to prevent a large stabilising 
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contribution from field-line bending, but -^^ is of order n. We also have -^Wip.Bi 



n. 



Because the system is axisymmetric, we need only consider a single Fourier mode in the 
toroidal angle, and take a ~ g-*"'/'. Then we have 







+ 



in 



0(1) 



(46) 



Later in Section |V] we will find that [iVV'.-Bil] = 0, so there will be zero contribution from 
the terms involving dJ^/dip. Rearranging the resulting equation leaves 



i d 

n dip 



Vip.Bi 



Bo.a + 



(47) 



Eq. |17|is our generalised criterion for marginal stability to the Peeling mode at high-n, with 
a given by Eq. 



E. Cylindrical limit 

Using Eq. Hi] (a = C,ipJo = JoVip.^), then for large aspect ratio (cylindrical) geometry 
a = Jo{RBp)^r and Eq. 1471 becomes 



RBp 



i dbr 
n dr 



Bq.Jq 
B. 7 



(48) 



Noting that K = B.V^r, taking ^ ~ gims-inc/,^ ^gj^g ± ^ j_ ^ with m ~ nq, this 



gives 

where we have also used that 
conditions (e.g. see Freidberg 



B„ 



m 



(49) 



or a cylinder q = . Finally the plasma vacuum boundary 



12|),of 



uq.B^ = Bq.V [no-U - [no.n UQ.nQ.V Bq 



(50) 

for a cylinder simplify to give = i?o-V^r = br- So we regain the usual condition for 
marginal stability to the Peeling mode in cylindrical geometry, of 



A^A, + J|| = 



(51) 



with A' 



r dbr 
br dr 



and A„=(l-^). 
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V. THE PLASMA- VACUUM BOUNDARY CONDITIONS 



The question now arises: what are the equivalent terms in the energy principle that 
correspond to our ordering for the Peeling mode? Before addressing this question, we first 
show that the plasma- vacuum boundary condition that is usually given in conjunction with 
the energy principle simply requires continuity of the normal component of the perturbed 
magnetic field, evaluated with the equilibrium normal to the surface at the equilibrium 
surface position. This is shown to agree with a simpler and more intuitive derivation that 
is given later. 

We start from the boundary condition usually given in conjunction with the energy 
principle (Freidbergjjj]), of 



= So.V ( no-e ) - ( ^0-^ ) ( ^0- (^o-V) fi, 



(52) 



with no the equilibrium unit normal. We will rewrite the expression in terms of Vip = noRBp, 
and then simplify the result. Firstly we multiply by RBp, and rearrange the equation to get 



V^p.BX = RBpBo.V (no.t) - (v^-f) (^-^V^.V^.VB, 
= Bo.V (RBpfio.tj - (no-tj Bo.V {RB^) 



V^.V [Bo.Vijj - Bo.Vij.VijVip 
5o.V - If 5o.V (RBp) - f« 



-Vip.Bo.V%l)Vip 



(53) 



-5n.V ''^ 



= Vtfj.Bi 

Alternately, V.-B = at the surface requires n.B = n.Bv, with n = no + ni, where no is 
the unit normal to the equilibrium magnetic flux surfaces (and hence also the equilibrium 
plasma surface) and ni is the unit normal perturbation from the equilibrium unit normal no 
(and hence also the perturbation from the equilibrium unit normal to the plasma surface). 
Then n.B = n.Bv requires 



no + ni) . 5o + 5i + eV no.Bo ] = {no + n{) .[B^ + B^ + no-^o^ 



(54) 



where ^ is the displacement of the plasma from its equilibrium position, and where all 
quantities are evaluated at their equilibrium positions. Similarly at equilibrium it is required 
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that 

no. Bo = uo-Bq (55) 

again with the quantities evaluated at their equihbrium positions. Therefore because 
uq-Bq = no-Bo = everywhere, then ^.V (^uo-Bq^ = ^.V (^uo-BX^ = 0, and Eq. [56] 
may be simphfied and rearranged as 

ni. (5o - + no. [b^ - = (56) 

Assuming there are no skin currents at equihbrium, then as shown previously in the main 
text, we have Bo = Bq , so Eq. [561 becomes 

no. (b^ - BX) = (57) 

or upon multiplying both sides by RBp 

V^. [Bi - BX) = (58) 

as above. 



VI. THE ENERGY PRINCIPLE 



Now we return to the relationship between our equations for marginal stability of the 
Peeling mode, and the energy principle. We start from the high mode number formulation 
for 5W = 6Wf + SWs + SWy, given in Connor et al[l^. Looking firstly at 6Ws, then later 



6Wv, we have 



SW^ = 71 (f) dx^J^Bkii 



R'Bll d 



J^B^ dip 



B.J 



B 



(59) 



with J^Bk^ = iJ^B.V. Note that Connor et al[lj] took ^ ~ g+«n0^ however here we have 
^ ~ g-«n</>^ which is reflected by n — > —n in the expression for 6W. We integrate by parts, 
noting that ^ ~ g-™'/' but ^* ~ e*"*^, and also replace J^Bk\\ with iJ^B.V to give 



SWs = -TT (b J^dx -B.V^ 



n 



R B i d 



Using B.Vil = V^J.Bl, then we get 



SWs = T^ <p Jxdx [ -^^-^1 



JyB'^ n dip 



R^Bli 1 d , 



B.J 

~B^ 



(60) 



B"^ n Jy dip 



BJ_ 



(61) 
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To obtain the vacuum solution for BY = W, we need to solve V^V = 0, which requires 



(62) 



1 and y ~ e we can rewrite this as 



{J^Vij.VV) + Lb. V + in 



1 / -> / \ n 

B.V + in— 1 l^- 



i?2 



i?2 



(63) 



Then using 5.V ~ 1, and taking the high-n limit, gives 



(64) 



which using = VV^, rearranges to give the result that for n ^ 1, 



B"^ dip 
The vacuum contribution to 6W is 

SWv = \ I dr 



J^Vip.Bl 



(65) 



bX 



(66) 



Using V'^V = 0, i?i = Vy, and Gauss' theorem, this may be written as an integral over the 
plasma surface, with 

5Wv = -\ I VdS.VV* (67) 



Using dS = ( ) (J^Bpdx) (Rdcj)), and integrating with respect to gives 



SWv = ~Tx <p JydxVVip.W* 



(68) 

Using the expression Eq. [65] for V at high-n, and the boundary condition Vijj.Bi = Vip.Bi 
(that Section |V] shows to be equivalent to the boundary condition that is usually given in 
formulations of the energy principle), we get 



SWv = f J^dx[Vi).Bi 
Inserting —l=i^ into SWv, we get 



R^Bl I I d f 



(69) 



SWs + 5Wv = r^§ J^dx (^) {vi,.Bl) [^^;^|, {j^Vip.B\ 



~B2 



(70) 
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Because V^.-Bf = VV^.-Bi, the terms involving dJ-^/dip will cancel. Then using the notation 
[I /I] to denote the difference between / evaluated just outside and just inside the plasma, 
gives 



5Ws + 



J^dx 



fi2 



i d 
n dip 



Vip.B] 



B.JC, 



(71) 



The term in {} is exactly Eq. |47| with a = Jq^^ as given by Eq. |4TJ Therefore marginal 
stability of the Peeling mode corresponds (at high-n) to taking the plasma's contribution to 
6W, 6Wf = 0, and then solving 6Ws + SWy = 0. This suggests we should define the high-n 
Peeling mode as a mode (represented by a trial function in this analysis), that allows us 
to neglect 6Wf compared to 6Ws + SWy (by for example being sufficiently localised), and 
whose subsequent stability is determined by 6Ws and 6Wv- 

VII. X-POINT PLASMAS 



Now we consider the stability of Peeling modes to the trial function considered by Laval 
et al 8|, that consists of a single Fourier mode with ^ = ^m(?/')e*'"^~*"'^, where ^ = ^ vdx' 
is the usual straight field-line poloidal coordinate. After taking Eq. [71], and using Vip.BY = 
Vip.Bi = B.V^^, we have 



SWs + 5W^ 



V 



7r 



■Jxdx - 



n 



52 



R^Bl 



Substituting the trial function into Eq. [72] gives 

|2 



6W 



where 



-2n- 



A 



2 \^rn\ 



Rq 



A AA' + J 



-]B.V^^ + B.J^^} (72) 



(73) 



m — nq 



J 



nq 

1 / „ IRo J.B 

2^f'^^WB;-w 



A' 



dlR^Bp 



dip 



(74) 
(75) 

(76) 



2tt J - ^K^B^ Vi).Bi 

and with dl = J^Bpdx an element of arc length in the poloidal cross-section, and Rq a typical 
measure of the major radius such as it's average for example. Note that because is a 
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Fourier component of ^.V?/' ~ ^{RBp), the dimensions of |^mP/-Ro are energy. Equation [731 
may easily be minimised for A (or equivalently, minimised with respect to choice of toroidal 
mode number), with 

A = (77) 
2A' ^ ^ 

giving 



If A is chosen to maximise the growth rate, then a similar but different value will be found. 
Similarly, there is no reason why there should not be a more unstable mode than the trial 
function we have considered. However, our primary interest is stability to the trial function 
that was found to be unstable by Laval et al Sj]. 

A calculation of 5W requires the evaluation of A' for a plasma equilibrium with a sep- 
aratrix. The calculation of A' is the main subject of the second part to this paper. As an 
introduction to this, we note that for a circular cross section plasma, an estimate for A' 
may be found by approximating the perturbation to the magnetic field near the edge of the 
plasma as being the same as for a vacuum, then solving Laplace's equation both inside and 
outside the plasma, and matching the solutions at the plasma-vacuum boundary. Then for 
a circular cross-section A' = —2m ~ —2nq. 

Observe that for J.V0 = 0, Bp. J = —I'Bp, for which 



2-K J R^Bp B^ 

and as discussed above we will have 6W — > as g — *• cxo. However, for J.'V<f) 7^ 0, then 
B.J = —Ip' — B^r and therefore 

^{Ro^)Hi (80) 
= ^0 (#) , 

For which if A' ~ —nq as is the case for a circular cross-section, then J^/A' would be of 
order —q, and 6W < 0, suggesting that the mode would be unstable. 

Although the sign of SW is usually taken to indicate whether a mode is unstable or not, 
the growth rate determines how unstable the mode is (i.e. how rapidly it develops). For 
example if our trial function = = ^'4'-^ had been .^m(^/')e*''"^ = RBp so 
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that it had dimensions of length as opposed to dimensions of length times RBp, then we 
would no longer have J q, despite our model only depending on the poloidal structure 
of the mode. The dependence of 5W on the normahsation of the plasma perturbation does 
not affect the calculation of the growth rate however, for which the consequences of the 
normalisation of the plasma perturbation will cancel. The growth rate is discussed later. 

Are there any reasons why a computer code might fail to find an unstable mode? One 
possibility is that the need for m ~ will require very high poloidal mode numbers as 
— > oo, and this could potentially prevent a numerical code from seeing the instability. Also 
important, is the need to consider the most unstable mode. Minimising 5W with respect 
to the toroidal mode number gave A = — J/2A'. If A' ~ —2nq and J ~ qRQ^J .B) / (for 
J.V0 7^ 0), this would require 



that is independent of q and the poloidal mode number. A final possibility is that A' 
might diverge more rapidly than J as we approach the separatrix. The second part to this 
paper calculates A' analytically, thereby avoiding the numerical problems associated with 
an X-point. 

Another, less obvious reason why an unstable mode might not be found in computer 
calculations is that despite 5W < indicating that it is energetically favourable for the 
mode to be unstable, the growth rate can still be vanishingly small. This possibility is 
explored next. 

VIII. THE GROWTH RATE 

So far we have only considered 5W , because its sign is usually presumed to be sufficient 
to indicate whether a mode is stable or not. The growth rate 7 for a mode with ^ ~ e'* is 
obtained from 7^ = —5W/\ J po\^\'^dr, where | / Po\^\'^dr is the kinetic energy term. Next 
we will estimate the kinetic energy term, so as to estimate the growth rate. The surprising 
result that we will find is that even if 6W < 0, indicating it is energetically favourable 
for an instability, the kinetic energy term can diverge so strongly that although it may be 
energetically favourable for a mode to be unstable, its growth rate is vanishingly small. An 
alternative complementary calculation to the one given below, with the same conclusions, is 




(81) 
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given in Appendix |X1 

To estimate J po\^\'^dr, we write 



^ R^Bl ^^^B"^ R^BIB"^ 



for which 



16 



10 



R^Bl 



52 



R^B^B^ 



(82) 



(83) 



It is convenient to write ^ in the form of Eq. [H21 so that we can use the results from a high-n 
ordering (e.g. see Webster and Wilson|19;]), that gives 



n n 



2 r>2 ^^tp 



and 



R^B^B^ 



B.VB^ 



(84) 



(85) 



Before continuing further we make some observations on the high-n ordering that for 
V.^ = usually leads to = -R^B"!^^. This analysis is used in the derivation of 6W 



p a-<p - n 

used in Section |VT] onwards, and to derive the equations solved by ELITE[9]. The ordering 

implicitly assumes that ^ ^ 7^75^5 which is the case in the plasma core where 7-75^ ~ 1, 



Jx dip 



Jx dip 



because ^ ~ n ^ 1. Whereas a sufficiently large n can always be found to ensure 
n7^^^ <C 1, and terms of this type will often be negligible order one contributions anyhow, 
future calculations would be improved by including them. For example would then 
become £1 = 1 -^p d(-h^i>) ^ the present we will continue to use the ordering employed 

an aip 
, that is also used to derive the equations solved by ELITE. 

Because the trial function that we consider consists of a single Fourier mode, Eq. [85] may 

be solved for ^b, with 



6 



; dip ) 



^ {im - inq) - 



givmg 



le 



B\ 



d_ 

dip 



{2p + 52 



+ 



9^ 

dip 



B.\/B^ 



B^ 



52 



(86) 



(87) 



with A 



^n2A2 + {b.WB^B^ 
{m — nq)/nq as before. Using Eq. [HUfor and substituting this and 1571 into Eq. 
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[83l gives 



\e = le 



i' 



+ 



dtp 



2 / B.WB^ 

B2 



+ 



I^Bp 2 A 2 



2B.VB^ 



Bl and ^ 



Noting that 

Bp 0, with A ~ 1/n as found previously, then we find 



Bp, whereas ( j^n^A'^ 



B.VB^ 



(88) 



B^/R^ as 



+ 



R^Bl 1 
52 n2 



(89) 



In other words, we may neglect the term in \^b\^ compared with the |^^|^ and terms. 
We will consider each of these terms in turn. Firstly 



(90) 



~ isj4,..j 



with ip = ips sX the plasma surface, and where we used[l8| 



?'(^) ^ 2^ / 



Bl 



J _ UOp 

dip 



dx 



~ 2-K J Bl 



(91) 



with (iZ = {J^Bp)dx, and at large aspect ratio V0.J is approximately a function of the 
poloidal magnetic fiux. Although we have not considered the radial structure of the mode 
in this paper, to estimate these terms we will adopt the ansatz that near the plasma's edge 
can be approximated by a power law, with 



im — io- 



(V^a-^)^ 



(92) 



where ifj = il^a Q-t the separatrix, and ifj = 4's the plasma surface. This is consistent 
with studies that do consider a mode's radial structure Q, [isl. We also use the result 



found here 



16j and elsewhere [20|], that for a conventional X-point (as opposed to the X- 



point produced by a "snowfiake" divertor[20[|). near a separatrix we have 



q ~ -go In 



(93) 
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for some constant go ~ 1- Under these assumptions 



q'M^di; ~ lev 



givmg 



Po 



\(: |2 



dr 



pgr 
B{B,) 



where we took IV (p J ~ B{Bp)/r, with r a measure of the plasma radius. 
Next we consider, 



Po- 



R^Bl 1 

Q2 ^2 



din 



dr 



52 



R Bpdl 



dip 



dip 



For the trial function = ^mi^P)^^'^^ with = - ^dx', we have 



din 



dip 



dip 



din 



dip 



+ \i. 



de 



dip 



I dip 



(94) 



(95) 



(96) 



(97) 



the last expression ignores terms in ^^|^ because either 



din 



or 



dip 



most strongly, so the largest contribution to the expression will be from either 
2 



will diverge 
2 

or 



\i. 



m 



2. Taking i^ as in Eq. [221 gives 





dim 


/ 


dip 



2/1 



dip ~ l^ol q 



Similarly, because ^ udx = 27rg and = - ^dx', we have 



de 



dip 



q' 1 f^du ^ , 
q q J dip 



that with m c::^ nq gives 



Jm'lir. 



2 


de 




dip 



n' 



'Iq"\i. 



Therefore, because ips ~ § B?Bpdl/ § dl we have 



Po 



R^Bl 1 
52 n2 



1p = i'a 



din 



d'4> 



(98) 



(99) 



(100) 



:ioi) 



Finally, taking V0.J 7^ and A' ~ —2nq (A' is calculated in the second part to this 
paper), then gives 



-6W 



7 



Jdrpolil'^ 



RqJ .B 



-ipsq' 



(102) 
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with 7a = B'^ I {pqR § dl) . Because g' — oo more rapidly than q as we approach the separa- 
trix, then for an outermost flux surface that is made increasingly close to that of a separatrix 
with ijjg — >• 'i/'a, we have that 7^ — and 

lnf,^)=-ilnf^l (103) 



with 7^1 the Alfven frequency, indicating that the growth rate 7 — as g' ^ 00. 
IX. SUMMARY 

This paper re-explores the stability of the Peeling mode for toroidal Tokamak geometry. 
It starts from a simple approach to Peeling mode stability at marginal stability in cylindrical 
geometry, then generalises this to toroidal Tokamak equilibrium. In the process of doing so 
we flnd a number of interesting results, namely 

1. At marginal stability, a plasma perturbation induces a skin current that is parallel 
and proportional to the equlibrium current at the edge, and proportional to the radial 
plasma displacement. 

2. For zero equilibrium skin current the usual plasma- vacuum boundary conditions 



(Freidberg[l2]), are identical to the requirement that uq.Bi = Hq.BY , with the quan- 



tities evaluated at the equilibrium position. 

3. The equilibrium conditions (force balance) for the Peeling mode at marginal stability 
and high toroidal mode number n, are identical to requiring 6Ws + SWy = 0, where 
6Ws and SWy are the surface and vacuum contributions to the energy principle's 
SW = 6Wp + 6Ws + SWy , with 6Wf the plasma's contribution to the energy principle 



(Freidberg 

This suggests the Peeling mode be deflned as a mode for which 6Wf <^ 6Ws + SWy ■ For the 
trial function used by Laval et al|8], that consisted of a single Fourier mode in straight fleld 
line co-ordinates, we flnd that the most unstable choice of A gives 6W = (^)^ i^^^- 
To evaluate SW for this model, it is necessary to know A' for a plasma cross-section with 
a separatrix and X-point at the plasma-vacuum boundary. Doing this without making the 
usual approximations (i.e. with a discretisation of space), as are usually made in numerical 
calculations, is the subject of the second part to this paperj 
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Finally we considered the growth rate, and found that even with SW < 0, the growth 
rate can be vanishingly small. This is because the kinetic energy term was found to diverge 
like q' as the outermost flux surface becomes increasingly close to a separatrix. When 
this divergence is sufficiently rapid (as would be the case for A' —2nq), then ln(7/7^) 
asymptotes to ln(7/7^) = —-^Inlipsg'/q), a result that may be compared with those from 
codes such as ELITE. 
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APPENDIX A: ALTERNATIVE DERIVATION FOR THE GROWTH RATE 

Here we provide an alternative derivation for the growth rate, to that given in Section 
IVIIII In the following we will. 



1. Continue to use the high-n ordering of Connor et al 
that 



14| . for which V.^ = requires 



2. We will use the arguments from Section IVlIIl that lead us to expect that d^^/dip ~ 
^l^ml ~ ^<l'\im\-, for the trial function of Laval et al[8,] with = ^m('^)e*™'^. 

3. As in Section IVlIII we will continue to assume that near the separatrix, 

(A2) 



for some constant go ~ !> with if) a the value of ip at the separatrix. We will also 
continue to assume that near the separatrix we can approximate l^ml"^)! as a power 
law, with 

where ipg < i'a is the value of tfj at the plasma-vacuum surface. 
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With the assumptions of[T],[2l and[3l we require that 

^-^ n P dtp 

^tR'B'^q'lU (A4) 

However, we must have ^± <^ 1 as ijj ^ ips, and therefore we require = ^o'^^^^^^^ with 
<^ 1 a constant, so that as ip ^ ipa we have ^± ~ ■Co ^ 1- Therefore we also have as a 
consequence that 

(A5) 

^0 as - V^a 



In the hmit where the plasma surface tends to a separatrix, with ipg ^ ipa, we then must 
have fj, ~ ,^0 ^^"""^^^ 0. Therefore 514^, for which fm is evaluated at = V's, has 

(A6) 



-> as 



Next we consider the growth rate. 

Using arguments from Section IVIIIt we expect 



I dr^^, (A7) 



Now with given by Eq. \AA\ we have 



fil lg±P _ r dld^d^R'^Bl |go|' 
J^^R^B^B^ J Bp B2 y^a-i^ J ^2 

P {B^) \ i'a 

where () denotes a poloidal average, and in the last line we used ipa ~ {Bp)R^. Hence using 
Eqs. |A6]and|A8]we find 

2 _ 5W 

/rfrpo||P 



^2 (^RoJ^\ ( 



pq 
n 
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with 7^ = {B"^) / pqRq § dl. Therefore as the outermost plasma surface more closely approxi- 
mates a separatrix with tps "ipa, we have that 7/7A 0. Note that because we have taken 



g ~ goln ( ^1 then q' ~ qo/i^a — '^), and hence (7/7a)^ ~ 1/?' as found Section IVIIII 
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